WDLEYKNOX  LIBRARY 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

GENERATION  AND  ANALYSIS  OF 
PARAMETERIZED  TERRAIN  FOR  LAND  COMBAT 

by 

Chung  Young  Kim 

March  19  77 

Thesis  Advisor:                 S.  H.  Parry 

Approved  for  public  release;  distribution  unlimited. 


T178062 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Otti  Enffd) 


DUDLEY  KNOX  UBMRV 
NAVAL  POSTGRADUATE  SCHOOL 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1.     KEF-OAT  NUMIER 

2.  GOVT  ACCESSION  NO. 

3.     RECIPIENT'S  CATALOG  NUMBER 

4.     TITLE  (and  Subtltlm) 

Generation   and  Analysis   of 
Parameterized  Terrain   for  Land  Combat 

5.     TYPE  OF  REPORT  *  PERIOD  COVERED 

Master's   Thesis; 
March    1977 

•  .    PERFORMING  ORG.  REPORT  NUMBER 

7.     AuTHOHf*; 

Chung  Young  Kim 

•      CONTRACT  OR  GRANT  NL-MBERfaJ 

t.     PERFORMING  ORGANIZATION   NAME   AND  ADDRESS 

Naval  Postgraduate   School 
Monterey,    California      93940 

10.    PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

1  1.     CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Postgraduate   School 
Monterey,    California      9  3940 

12.     REPORT  DATE 

March   19  77 

13.     NUMBER  OF  PAGES 
59 

14.     MONITORING   AGENCY  NAME  ft    ADDRESSf//  dlllaranl  from  Controlling  Ol'llea) 

18.    SECURITY  CLASS,  (ol  thla  rdport) 

Unclassified 

ISa.     DECLASSIFI  CATION/ DOWN  GRADING 
SCHEDULE 

16.     DISTRIBUTION   STATEMENT  (ol  tr\li  Rmpori) 

Approved   for  public   release;    distribution  unlimited. 

17.     DISTRIBUTION  STATEMENT  (ol  tfta  aaatract  «ii»r««  In  Block  20,  II  dlllaranl  from  Raport) 

!•.     SUPPLEMENTARY  NOTES 

Parametric  Terrain 

20.     ABSTRACT  (Contlmta  on  rwwM  aid*  it  naeaaaarr  «•»»*  Idmntlly  by  Aloe*  m— BwJ 

Many   current  army   combat   simulation  models    require 
digitized   terrain   inputs    for  their  execution,    resulting 
in   a  costly   and  time-consuming  operation.      A  methodology 
to   generate   parametric   terrain   as   a   function   of   the   number 
of  hill  masses   and  the  height,    location,    spread,    and  slope 
of  the  hills   is   presented  in   this   paper.      Modified  Bivariate 
Normal    (MBVN)    distribution   function   and   a   linear   function 

DD  ,* 

(Page    1) 


Fa-:M7,  1473 


EDITION  OF  1  NOV  61  IS  OBSOLETE 
S/N   0102-014*6601  I 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Dmta  tntarad) 


UNCLASSIFIED 


ffnOJ'WTV   CLASSIFICATION   OF    TmS  PtflErt^w   n»f«  Enfr»J 


(20.   ABSTRACT   Continued) 

cutting  the  MBVN  function  are  employed  to  achieve  the  terrain 
representation.   Analysis  of  various  terrain  parameter  values 
was  conducted  and  the  significant  trends  resulting  from  the 
study  are  described.   Finally,  recommendations  for  future 
enhancement  of  the  model  are  given. 


DD     Form       1473  -— ,«-  ■*  __._„_„_, 

1  Jan  73  UNCLASSIFIED 

S/N    0102-014-6601  0  security  classification  of  this  PAcer****  o«r«  ««»•'•<*) 


Approved  for  public  release;  distribution  unlimited 

Generation  and  Analysis  of 
Parameterized  Terrain  for  Land  Combat 

by 

Chung  Young  Kim 
Lieutenant  Colonel,  Korean  Army 
B.S.,  ROK  Military  Academy,  1963 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
March  1977 


ABSTRACT 

Many  current  army  combat  simulation  models  require 
digitized  terrain  inputs  for  their  execution,  resulting 
in  a  costly  and  time-consuming  operation.   A  methodology 
to  generate  parametric  terrain  as  a  function  of  the 
number  of  hill  masses  and  the  height,  location,  spread, 
and  slope  of  the  hills  is  presented  in  this  paper.   Modi- 
fied bivariate  normal  (MBVN)  distribution  function  and  a 
linear  function  cutting  the  MBVN  function  are  employed  to 
achieve  the  terrain  representation.   Analysis  of  various 
terrain  parameter  values  was  conducted  and  the  significant 
trends  resulting  from  the  study  are  described.   Finally, 
recommendations  for  future  enhancement  of  the  model  are 
given. 
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I.   INTRODUCTION 

Many  current  army  combat  simulation  models  require 
digitized  terrain  inputs  for  their  execution.   This  process 
is  very  costly  and  requires  a  significant  amount  of  computer 
storage.   In  addition,  computation  of  line-of -sight  and 
mobility  require  a  great  deal  of  computer  running  time. 

The  current  methodology  utilized  to  represent  terrain 
is  also  questionable  from  a  statistical  point  of  view. 
A  Vector  Research  Incorporated  (VRI)  study  on  terrain  line 
of  sight,  concluded  that  "...  present  and  past  army  study 
results,  based  on  the  analysis  of  combat  results  on  a  very 
limited  sample  of  terrains,  may  have  been  determined  by 
the  terrain  selection  process  and  not  by  the  actual  weapon 
system  or  force  design  differences."  More  specifically 
it  stated  the  following: 

A.  "There  is  extreme  sensitivity  in  combat  model 
results  as  the  scenarios  (terrain  and  movement  assumption) 
are  varied,  even  when  variation  is  within  a  class  of 
scenarios  chosen  for  their  a  priori  equivalence." 

B.  "This  sensitivity  can  be  slightly  reduced,  but 
remains  extreme  (with  probabilities  of  win  estimable  only 

within  plus  or  minus  25%)  even  when  battle  results  are 

«1 
used  to  redesign  scenarios." 


Farrell,  Robert  L.  and  Richard  J.  Freedman,  Investigation 
of  the  Variation  of  Combat  Model  Predictions  with  Terrain  Line 
of  Sight,  Vector  Research  Incorporated,  1975. 


These  results  imply  that  sufficient  replications  of  each 
type  of  terrain  should  be  run  in  order  to  reach  a  satisfactory 
statistical  level  of  significance.   For  VRI's  analysis  at 
least  fifty  replications  of  each  type  of  terrain  were  used. 

In  view  of  the  problems  stated  above,  research  was  under- 
taken  by  Major  Christopher  Needels   to  develop  a  new  method^ 
ology  for  the  representation  of  terrain.   In  particular,  a 
modified  bivariate  normal  distribution  function  was  utilized 
to  generate  hill  masses  in  a  wide  variety  of  configurations. 
Line  of  sight  and  movement  routines  were  also  developed  to 
evaluate  various  intervisibility  parameters  over  many  terrain 
configurations  for  very  low  cost. 

This  thesis  represents  an  extension  of  the  models 
developed  by  Needels  to  generalize  the  terrain  generation 
routines.   In  particular,  methodology  was  developed  to 
generalize  the  hill  mass  generation  from  symmetric  hill 
shapes  to  non-symmetric  configurations.   The  methodology 
developed  for  this  generalization  is  described  in  Chapter  III. 

In  order  to  evaluate  the  effects  of  the  number  of  hills, 
hill  height,  and  slope  on  various  intervisibility  parameters, 
an  experiment  was  designed  to  measure  these  effects.   A 
description  of  the  design  methodology  and  results  of  the 
analysis  are  given  in  Chapter  IV.   Finally,  the  conclusions 
drawn  from  the  current  analysis  and  recommendations  for 
future  research  are  presented  in  Chapter  V. 


2Needels,  Christopher  J.,  Parameterization  of  Terrain  in 
Army  Combat  Analysis,  Master's  Thesis,  Naval  Postgraduate 
School,  Monterey,  California,  1976. 


II.   REPRESENTATION  OF  TERRAIN 

A.   HIGH  RESOLUTION  COMBAT  MODEL  REPRESENTATION 

The  classic  approach  to  representing  terrain  in  land 
combat  models  is  to  utilize  digitized  terrain  evaluation 
data  from  Waterway  Experiment  Station  (WES)  for  the 
particular  area  of  interest. 

The  Dynamic  Tactical  Simulator  (DYNTACS)  model  serves 
as  a  good  example  of  the  state  of  the  art  in  terrain 
representation  and  use.   DYNTACS  is  a  two-sided  dynamic 
model  of  battalion  level  combat  representing  the  detailed 
interaction  of  elements  in  a  combined  arms  environment. 

DYNTACS  utilizes  100  meter  grid  squares  for  which 
macro-terrain  elevations  are  determined  from  WES  tapes. 
The  simulation  divides  each  grid  square  diagonally,  thus 
providing  a  series  of  adjoining  triangular  terrain.   The 
entire  battlefield,  therefore,  is  represented  as  a  surface 
of  equal  size  triangles  which  vary  in  slope  depending  on 
the  elevation  at  their  corners. 

Needels  (pp.  9-10)  provides  a  concise  description  of 

the  DYNTACS  methodology  for  line  of  sight  calculation  as 

follows. 

From  the  macro-terrain  data,  line  of  sight 
between  any  two  opposing  elements  is  computed. 
This  is  accomplished  by  first  computing  the 
angle  between  the  horizontal  and  a  straight 
line  drawn  between  an  observer  and  target 
(O-T  Line) .   The  program  then  conducts  a 
search  of  the  terrain  along  the  path  of  the 
O-T  line  to  see  if  any  macro-terrain  is 


higher  than  the  0-T  line  itself.   This  is 
accomplished  by  comparing  the  angle  of  the 
0-T  line  with  the  angle  above  horizontal 
of  the  Observer-Terrain  line.   If  the  latter 
angle  is  larger,  there  is  no  intervisibility . 
Over  the  duration  of  a  battle  with  numerous 
elements  this  calculation  may  be  made 
thousands  of  times.   Consequently,  not  only 
is  the  time  to  prepare  the  tapes  high,  but 
also  the  time  to  compute  lines  of  sight  once 
the  terrain  is  input  to  the  model. 
Considerable  effort  by  developers  and  users 
of  this  model  has  been  expended  to  stream- 
line this  subroutine. 


B.   SYMMETRIC  REPRESENTATION 

3 

The  basic  models  developed  by  Needels   (hereafter  referred 

to  as  SIMTER)  which  provide  the  basic  foundation  for  the 
methodology  developed  in  this  thesis  are  described  in  this 
section.   The  basic  motivation  for  SIMTER  was  to  develop  a 
random  terrain  representation  which  could  be  replicated 
quickly  but  at  the  same  time  be  representative  of  a  particular 
"type"  of  terrain.   The  approach  adopted  for  SIMTER  generated 
a  variable  number  of  hill  masses.   The  shape  and  height  of 
the  hills  could  be  varied,  but  each  hill  was  symmetric,  since 
it  was  generated  directly  from  the  modified  bivariate  normal 
(MBVN)  density  function. 

SIMTER  has  two  parts.   The  first  part  is  the  main 
program  which  creates  terrain  and,  at  the  discretion  of  the 
user,  plots  both  a  three-dimensional  drawing  and  a  contour 
map.   Used  strictly  for  terrain  generation,  it  can  be  used 
as  a  preprocessor  for  other  combat  models  by  producing  grid 
points  and  elevations  similar  to  those  provided  on  computer 


3Needels ,  Parameterization  of  Terrain  in  Army  Combat 
Analysis . 


tapes  by  WES.  The  second  part  of  SIMTER  is  a  movement/line 
of  sight  routine  which  moves  a  target  along  specified  routes 
across  the  generated  terrain. 

C.   A  NEW  REPRESENTATION  METHODOLOGY 

The  methodology  developed  in  this  thesis  provides  for 
the  generation  of  hill  masses  utilizing  MBVN  density 
function  to  represent  masses  that  are  not  symmetric.   This 
development  provides  a  significant  increase  in  flexibility 
to  represent  hill  masses  more  realistically.   The  mathematical 
development  of  this  methodology  is  presented  in  the  next 
chapter,  after  which  actual  results  from  executing  the  model 
are  given. 
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III.   THE  EXTENDED  MODEL 

A.   UNSYMMETRIC  TERRAIN 

The  unsymmetric  slope  of  terrain  can  be  generated  by 
a  cutting  method  using  the  MBVN  as  follows: 


_i_  [(_x,2  _  2  (_5)<_JL)  4  (-^-Y-)2] 
A  ( 1-  p  )      x  x     y       y 


where    MH:  maximum  height  of  hill 

Z  :  height  of  hill 

]i    rV  •    spread  of  hill 

p  :  ellipse  factor  of  hill 


The  unsymmetric  slope  of  the  hill  can  be  represented  by  the 
vertical  value  of  the  MBVN  subtracted  from  the  vertical  value 
of  a  linear  function  as  shown  in  Figure  1  by  the  value  HI. 


/ 

\ 

A  \ 

/hi  \ 


FIGURE  1.   UNSYMMETRIC  SLOPE  GENERATION 
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Two  basic  questions  must  be  addressed: 

1.  Where  do  we  start  the  cutting  of  the  MBVN  function? 

2.  How  do  we  compute  the  initial  cutting  position? 
If  M  is  the  direction  of  the  cutting  line  (Line  A) ,  then 
the  cutting  plane  is  through  the  location  of  the  center  of 
the  hill  (Ux/Py)  as  shown  in  Figure  2. 


FIGURE  2.   CUTTING  LINE  DETERMINATION 

If  the  initial  point  of  cutting  is  (XIN,YIN)  and  HI  is 
height  at  the  initial  point,  this  point  is  on  the  cutting 
plane  and  also  on  the  MBVN  distribution  function.   Thus,  if 
we  consider  the  graph  in  two  dimensions  as  shown  in  Figure  2, 
Line  B  is  tangent  to  the  circle  at  the  point  that  the  vertical 
value  on  the  MBVN  function  is  HI  and  the  slope  of  Line  B 
is  -1/M.   If  the  initial  point  of  the  cutting  is  (XIN,YIN) 
this  point  is  on  Line  B,  Line  A  and  the  circle  at  a  height 
HI  on  the  MBVN  function.   We  can  compute  the  initial  point 
(XIN,YIN)  as  follows. 
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1.   p=0  Case 


1   XIN^X  2    YIN^V  2 
HI   =   MH*EXP}-  •=■[( -)*    +  ( £)Z] 


(1) 


-2  ln(HI/MH)   = 


XIN-u        YIN-y 
(— V  +  (— Z)z 


(2) 


From  (2)  ,  the  derivative  is  computed  as  follows 

XIN-y 

( 

dYIN 


( — A 


dXIN 


YIN-y 

( 2~) 


1 

M 


(3) 


From  (2)  and  (3) 


XIN  = 


y   ±  a 
Mx    x 


-2  ln(HI/MH) 
l+M2(ay/ax)2 


(4) 


YIN   = 


y   ±  a  M 
My     x 


-2  ln(HI/MH) 
/ay) 2  +  M2 


(5) 


If  ox  =  oy    ,    then 


XIN  = 


(6) 


.  -2  In  HI/MH) 
YIN      yv  ±  a  M  2 

y   /   i  +  w 


(7) 
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If  the  hill  is  cut  from  left  to  right  and  a     =  a    ,    then 


XIN  = 


-2  ln(HI/MH) 
X     X  '    1  +  M2 


(8) 


YIN   = 


_  a   M/-2  ln(HI/MH) 

y        y   '   1  +  M2 


(9) 


If  the  hill  is  cut  from  right  to  left  and  a   =  a  ,  then 


XIN   = 


y   +  a 
x    x 


-2  ln(HI/MH) 


1  +  M' 


(10) 


YIN   = 


y  +  a   M 

y   y 


-2  ln(HI/MH) 


1  +  M' 


(11) 


Thus,  if  the  direction  of  cutting  and  HI  are  selected,  then 
the  starting  point  can  be  computed.   Care  must  be  taken  in 
the  selection  of  the  value  of  HI.   If  HI  is  too  small,  then 
the  hill  will  have  an  unrealistic  shape.   The  value  of  HI 
affects  the  shape  of  the  hill  as  shown  in  Figure  3. 


FIGURE  3.   EFFECT  OF  HI  ON  THE  HILL  SHAPE 
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2.   0  <  p  <  1  Case 


,      XIN-y v  0      XIN-y    YIN-y 

HI   =   MH  *  EXP{ ^-o-[  ( -)   -  2  ( -)  ( £) 

2(l-p  )     x  x      y 


YIN-y   0  N 

+  (-^'i  (12> 


9               XIN-y   9       XIN-y    YIN-y 
-2(1  -p^lnCHI/MH)  =  (— *-r  -  2  (— -)  (— *) 


a  a       a 

x  x       y 


YIN-y 
+  (      Y)2         (13) 

y 


(14) 


From  (13)  the  derivative  is  computed  as  follows 


,    XIN-y        YIN-y 

-  :r-  [  (  „     x)  +  (  =  y)  1 

ax     ax  Gy 

dYIN      1 

dXIN  "  "  M 

XIN-y       YIN-y 


From  (14) 


XIN-yv       a   +  Mpa    YIN-y 

( S)   =   -* £  (      v)  (15) 

ax        May  +  pax     y 
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Let   C  = 


a     +  Mpa 
x      x 

Ma   +  pa 
y     x 


From  (13) ,  (15)  and  (16) , 


-2(1 -p*)  ln(HI/MH)   =   ( 


~     ,   XIN-y   0 
C      x 


(17) 


Let   D  = 


1  -  2fi.+  4r 

C    C2 


(18) 


Then 


XIN  =   u   ±  a 
Kx    x 


-2(1  -  p*)  ln(HI/MH) 


(19) 


Also,  let 


E   =   C   +  2pC  +  1 


(20) 


Then, 


YIN   = 


yy  *  ay 


-2(1  -q)    ln(HI/MH) 


(21) 
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FIGURE  4.   CUTTING  DIRECTION 

If  the  cutting  direction  goes  through  the  longer 
axis  as  shown  in  Figure  4 ,  then 


r,  2    „  2 

cot  2a   =  -A, %— 

-2pa  a 
x  y 


(22) 


The  cutting  direction,  a,  is  determined  by 


»      2     rr      2 

,        a   -  a 

a  =  T  arccot  ?     * 
2         -2pa  a 
x  y 


(23) 


If  the  cutting  direction  goes  through  the  shorter 
axis  as  shown  in  Figure  4,  then 


6   =   90°  -  a 


(24) 


18 


Next,  the  height  of  intersecting  hill  masses  must 
be  determined  as  illustrated  in  Figure  5. 


FIGURE  5.   INTERSECTING  HILL  MASSES 

HB   defines  height  of  hill  "B"  after  cutting. 
HA  defines  height  of  hill  "A"  after  cutting. 

Height  of  hill  "A"  is  greater  than  that  of  hill  "B" 
at  the  point  XA  but  height  of  hill  "B"  is  greater  than  that 
of  hill  "A"  after  cutting.   If  there  is  more  than  one  hill 
located  close  to  each  other,  those  hills  must  be  adjusted 
by  comparing  the  height  after  cutting  each  one .   The 
computation  of  the  height  of  the  terrain  at  any  point  is 
described  in  the  following  section. 

B.   TERRAIN  HEIGHT 

The  height  of  the  hill  which  is  generated  by  the  MBVN 
distribution  corresponds  to  the  height  of  a  hill  which  is 
generated  by  the  cutting  method.   The  height  of  a  hill 
at  any  point  is  computed  as  follows: 
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The  height  of  the  MBVN  density  function  is  given  by 


x-y 


1       •*•  M„  o       x-y   y-u      v-u 
Z   -   MH*EXP  {-  — 1     [(-^S)2  -  2p(_£)(L-*,  +(Li)2l> 


2(l-p*) 


a   "a 
x    y 


The  height  of  the  linear  function  is 


(25) 


ZCUT   =   (X-XIN)tan6+  (Y  -  YIN)  tan  <J> 


(26) 


where     9  =  cutting  angle  of  the  x-axis 
<$>   =  cutting  angle  of  the  y-axis 

The  height  of  the  hill  after  cutting  is 


ZH   = 


f    ° 

0 

Z  -  HI 

Z  -  ZCUT  -HI 
0 


if   Z  <  HI   and  ZCUT  <  0 

if   Z  <  HI   and  ZCUT  >  0 

if   Z  >  HI   and  ZCUT  >  0 

if   Z  >  HI   and  ZCUT  >  0 

and  Z  >  ZCUT+HI 

if   Z  >  HI   and  ZCUT  >  0 

and  Z  <  ZCUT+HI 


If  more  than  one  hill  exists  at  a  point  on  the  terrain,  the 
greatest  height  after  cutting  is  used. 

C.   MOVEMENT  LINE  OF  SIGHT 

The  procedure  of  computing  the  line-of-sight  between 
any  two  points  on  the  terrain  is  described  in  this  section. 
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If  the  position  of  the  observer  is  (OX,OY)  and  the  position 
of  the  target  is  (X,Y) ,  then  the  slope  of  the  Observer- 
Target  (0-T)  Line  is  given  by 


ML   = 


X  -  OX 
Y  -  OY 


(28) 


If  the  intersection  point  of  a  hill  and  the  0-T  Line  is 
(XI, YI) ,  then  HI  and  the  0-T  Line  are  perpendicular  as 
shown  in  Figure  6 . 


\ 


CXL.yz) 


hA*(.MtjMy)   . 


/   / 


cx,y) 


V  V 


S 


Thus  , 


FIGURE  6.   INTERSECTION  OF  0-T  LINE 
WITH  HILL  MASS 


XI   = 


ML   •  OX  -  ML  •  OY  +  y  +  ML  .  y 


1  +  ML 


(29) 


YI   =   ML (XI  -  OX)  +  OY 


(30) 
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If  there  is  more  than  one  hill  on  the  0-T  Line,  the 
dominant  hill  (after  cutting)  is  used  to  compute  the  percent 
of  the  target  covered  relative  to  the  observer.   Target 
movement  along  the  surface  slope  of  the  hill  is  computed 
as  follows: 

,       x-y   0      x-y   y-u 

ZH   =   MH*EXP  { ±— 5-  [(-~)2    ~    2p(— -£)(-— *) 

2(lV)     Qx  ax    ay 

y 

-  (X-XIN)  tan  0  -  (Y-YIN)  tan  <|>         (31) 


3ZH       -MH    *~%    Pty-V 


x-y   9     x-y   y-y      y-u  2 
*EXP{   -1  ,  [(-tt^)2  '2q(—^)(—^)    +  (-^)2]>  -  tan  0 
2(l-p^)    x        x    y       y 

(32) 


x-u  p (y~u  )  • 

3ZH      -MH   r,__x,  __X-1 

"8v~~  =  7  [( 2")  "  a  a    j 

d*      l-o     a  x  y 


*EXP{-Z±^.[{-^JL)2   -2p  (—£)(—*)  +  (-g-*)  21>-  tan  ♦ 


2(l-p2)   ax       ~x    -y       y 


(33) 
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If  DEL  is  the  direction  of  target  movement,  then 


DEL  =  |P  cos  r  +  |2£  sin  r  (34) 


If  the  speed  of  the  target  is  constant  "V"  on  the  surface, 
then 


velocity  =  ■  (35) 


/  1  +  (DEL)2 


If  the  speed  of  the  target  is  dependent  on  the  slope  of 
the  surface,  then 

velocity   =   V  +  W(DEL)  {36) 

where  W  is  the  speed  factor  of  slope. 
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IV.   TEST  AND  EVALUATION  OF  MODEL  RESULTS 

A.   DESIGN  AND  CUTTING  ANGLE  METHODOLOGY 

The  methodology  for  parametrically  describing  terrain 
as  a  function  of  the  number  of  hills,  the  center  and  spread 
of  each  hill,  and  the  slope  was  developed  in  the  previous 
chapters.   In  this  chapter  the  effect  of  the  number  of  hills, 
height  of  hills,  and  the  slope  on  various  intervisibility 
measures  is  investigated.   Three  random  terrains  were 
generated  (with  each  terrain  having  six,  twelve,  and  eighteen 
hills)  resulting  in  nine  basic  terrain  configurations. 
For  each  configuration  three  average  hill  heights  (100, 
200,  and  300  meters)  were  utilized. 

Finally,  for  each  of  the  27  resulting  configurations, 
four  slope  cutting  methods  were  used  as  follows: 

1.  Down  slope-cut  from  observer  side  to  target  side. 

2.  Side  slope-cut  from  right  to  left  as  viewed  from 
the  target. 

3.  Up  slope  —  cut  from  target  side  to  observer  side. 

4.  Both  slope-symmetric  slope  as  viewed  from  any 
position. 

A  total  of  108  unique  terrain  configurations  were 
generated  for  the  investigation.   For  each  of  the  config- 
urations, the  mean  and  standard  deviation  describing  the 
center  and  spread  of  each  hill  is  specified.   In  addition, 
the  cutting  angle  for  determining  slope  must  also  be 
specified. 
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If  the  cutting  angle  is  constant,  then  the  hill  can  be 
described  as  follows: 

HILL  =   F(Peak,  Position,  Standard  Deviation,  Slope) 

In  the  constant  cutting  angle  case,  only  the  slope  is 
variable,  with  peak,  standard  deviation,  and  position 
constant.   In  other  words,  if  the  slope  varies  in  Figure  7, 
the  standard  deviation  of  the  resulting  distribution  also 
varies. 


A  LINE- 


FIGURE  7.   CUTTING  ANGLE  AND  SLOPE 

The  objective  of  the  methodology  is  to  maintain  a  constant 
standard  deviation  of  the  distribuition  function   resulting 
from  each  cutting  angle,  so  that  only  the  slope  of  the  hill 
varies.   Otherwise  the  determination  of  how  slope  affects 
the  intervisibility  parameters  would  not  be  possible,  since 
it  would  be  confound  with  the  standard  deviation  (spread) 
of  the  hills. 
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The  problem,  therefore,  was  to  determine  a  method  to 
use  various  cutting  angles  and  maintain  a  constant  standard 
deviation.   After  extensive  investigations  into  various 
approaches  to  this  problem,  no  method  was  found  which  would 
mathematically  guarantee  a  constant  standard  deviation. 
The  method  shown  in  Figure  8  was  adapted  as  the  best 
approximation  method  available. 

\\tv£    ft 


L/fO^r    C 


FIGURE  8.   CUTTING  ANGLE  METHODOLOGY 

In  Figure  8,  HT  is  the  maximum  height  of  the  hill  after 
cutting  with  Line  A.   Note  that  the  cutting  angle  is  now 
specified  in  reference  to  the  center  line  of  the  original 
bivariate  normal  hill  mass  distribution.   In  order  to 
generate  a  symmetric  hill  mass  with  approximately  the  same 
standard  deviation  (Line  C)  the  cut  is  made  maintaining 
a  constant  maximum  height,  HT,  as  in  Figure  8. 

Several  test  runs  made  using  this  methodology  indicated 
that  the  resulting  standard  deviations  were  very  nearly 
equal  for  cutting  angles.   However,  very  little  change  in 
slope  results  for  cutting  angles  greater  than  twenty  degrees 
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B.  OBSERVER- TARGET  POSITIONS 

For  each  terrain  an  observer  position  and  initial  target 
position  were  selected.   Two  routes  from  the  target  to  the 
observer  position  were  selected  in  accordance  with  accepted 
tactics  of  advance  routes.   Selected  terrains  indicating 
the  observer-target  position  and  selected  routes  are  shown 
in  Figures  12  through  17.   in  order  to  assure  that  the 
"number  of  hills"  factor  was  properly  analyzed,  the  same 
routes  of  approach  on  each  generated  terrain  were  used, 
independent  of  the  number  of  hills. 

C.  TARGET  MOVEMENT 

Two  cases  of  target  movement  were  considered  in  the 
analysis . 

1.  Constant  Speed 

The  target  moves  along  the  route  at  a  constant  speed 
of  six  meters  per  second,  independent  of  hill  slopes.   The 
constant  speed  runs  serve  as  the  base  case  for  measurement 
of  the  percent  of  time  intervisible  along  the  route. 

2.  Variable  Speed 

The  speed  of  the  target  is  a  function  of  the  terrain 
slope  being  traversed.  The  speed  of  the  vehicle  is  computed 
by  the  relationship  as  follows: 

velocity  =   6.0  +  4.0 (slope) 

This  slope  is  the  directional  derivative  value  at  the 
instantaneous  target  position  on  the  hill. 
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D.  TARGET  INTERVISIBILITY 

Five  dependent  variables  describing  target  intervisibility 
are  computed  for  each  advance  route  as  follows: 

1.  Percent  of  time  intervisible  along  the  route, 

2.  The  number  of  intervisibility  segments  along  the  route, 

3.  The  average  intervisibility  segment  length, 

4.  The  average  intervisibility  segment  time,  and 

5.  The  total  length  of  intervisibility  over  the  route. 
Intervisibility  parameters  are  computed  by  utilizing  the 
line  of  sight  algorithm  at  fixed  time  intervals,  TL,  during 
the  movement  run.   For  this  analysis,  TL  was  set  equal  to 
one  second.   A  sample  output  for  one  compute  movement  trace 
is  given  in  Figure  9. 

E.  EXPERIMENTAL  DESIGN 

The  objective  of  the  analysis  is  to  determine  the  effects 
of  the  number  of  hills  (A) ,  the  height  of  hills  (B)  and  the 
slope  of  the  hills  (R)  on  target  intervisibility.   To 
accomplish  this  objective,  an  experiment  was  designed  as 
described  by  the  following  design  matrix. 

SLOPE 
#  of  Hills   Height     Down    Side     Up     Both 

100 

6        200 

300 

100 

12         200 

300 

100 

18        200 

300 
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FIGURE  9.   OUTPUT  OF  DATA 
Definition  of  Variables 


N 

XMEL,YMEL 

HITE 


Number  of  hills 
Spread  of  hills 
Maximum  height  of  hill 


RUFX,RUFY,RUFP  Different  amount  of  the  spread  and  height 

among  hills 


INCX,INCY 
XMAP,YMAP 
XGRID,YGRID 
RIGHT 

XSTD,YSTD 

XSLOP,YSLOP 

DIRECTION 

XINIT,YINIT 

HEIGHT 

RX,RY 

TIME 


Grid  increment  for  the  terrain  matrix 

Size  of  the  grid  system  in  meters 

Center  position  of  hill 

Maximum  height  of  hill  in  modified  bivariate 
normal  distribution 

Spread  of  hill  (standard  deviation) 

Cutting  angle 

Cutting  direction 

Initial  position  of  cutting  plane 

Hieght  of  hill  at  the  center  position  on  MBVN 
after  cutting 

Coordinates  of  approaching  route  of  target 

Movement  time 


MAP  COORDINATES   Target  position  at  given  time 


ELEVATION 
SPEED 
VISIBLE  ■ 
HLIN 


Target  elevation  on  the  terrain 

Velocity  on  the  terrain  at  the  given  time 

1  means  visible  and  0  means  invisible 

The  different  height  of  0-T  line  at  the 
maximum  height  of  surface  of  the  terrain  on 
the  0-T  line,  and  the  lower  height  between 
observer  and  target 


DZ  Maximum  height  of  surface  of  the  terrain  on  the 

0-T  line 

KL  The  hill  which  is  maximum  height  between  the 

observer  and  target 

DLIN        Difference  maximum  height  of  surface  between 
observer  and  target,  and  lower  height  between 
observer  and  target. 
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The  appropriate  linear  model  is  as  follows: 


Yijk    =  V   +  A.  +  B.  +  Rk  +  AB.  .  4-  AR.k  +  BR.k  +  ABR 

:ijk£ 


J     K     ID      ik      jk       13k 

+  e 


where 


Yijk£  =  intervisit>ility  dependent  variable  measure 
y     =  grand  mean  of  the  population 

A     =  number  of  hills,  where 

1-6  hills 
i  =   2-12  hills 
3-18  hills 

B     =  average  height  of  hills,  where 

1-100  meters 

j  =   2-200  meters 

3-300  meters 

R     =  direction  of  slope  of  hills,  where 

1-down  slope 
.     2-side  slope 
~  3-up  slope 

4-both  (symmetric) 

Within  each  cell  of  the  matrix,  there  are  six  replications 

consisting  of  two  routes  on  each  of  the  three  basic  terrain 

configurations . 

For  the  first  treatment  levels  of  slope  defined  above, 

a  cutting  angle  of  fifteen  degrees  was  used,  the  symmetric 

case  was  generated  with  a  cutting  angle  of  zero  degrees 

using  the  methodology  previously  described.   Four  analyses 

were  conducted  as  representative  of  the  large  number  of 

possible  sitatuations .   The  measurements  selected  were 

1.  Percent  of  time  intervisible  for  constant  speed 
and  for  variable  speed,  and 

2.  Number  of  intervisibility  segments  for  constant  speed 
and  for  variable  speed. 
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The  results  for  each  selected  analysis  are  described  in 
the  following  section. 

F.   ANALYSIS  OF  RESULTS 

!•   Percentage  of  Intervisibility  Time 
a.   Constant  Speed 


/  1/8  m 

(x    Transgeneration  Data) 


SUMMARY  TABLE  -  ANOVA 

Source        DF  SS        MS        F  Ratio 

#  of  hills        2  0.0437  0.0218       23.6417 
A 

height  of  hill    2  0.0117  0.0058       6.3209 
B 

slope  of  hill     3  0.0387  0.0129      13.962 
R 

AB  interaction    4  0.0006  0.0081       0.1527 

AR  interaction    6  0.0026  0.0004       0.4759 

BR  interaction   6  0.0016  0.00027      0.2943 

ABR  interaction  12  0.0012  0.0001       0.1114 

Error          179  0.1655  0.0009 


Each  of  the  three  main  effects  above  are  significant 
at  the  0.01  level,  but  no  interaction  terms  are  significant. 

The  following  table  gives  the  mean  and  standard  deviation 
for  each  treatment  level  of  slope. 
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Slope 

Mean 

Std.  Dev. 

down 

0.6130 

0.1539 

side 

0.5214 

0.1701 

up 

0.6926 

0.1629 

both 

0.5573 

0.1595 

Note  that  the  up  slope  has  the  greatest  percentage  of  inter- 
visibility  time;  down  slope,  both, and  side  slope  following 
in  that  order. 

The  following  table  gives  the  mean  and  standard 
deviation  for  the  treatment  levels  of  number  of  hills  and 
height  of  hills  for  the  constant  speed  case. 

#  of  Hills Mean Standard  Deviation 

6  0.6976  0.1620 

12  0.5462  0.1620 

18  0.5272  0.1490 

Height  of  Hill Mean Standard  Deviation 

100  0.6475  0.1775 

200  0.5844  0.1684 

300  0.5556  0.1628 
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b.   Speed  which  Depends  on  Slope 


t    3/8  „ 

(x    Trans generation  Data) 


SUMMARY  TABLE  -  ANOVA 


Source 

DF 

SS 

MS 

F  Ratio 

#  of  hills 
A 

2 

0.5152 

0.2576 

34.286 

Height  of  hill 
B 

2 

0.2041 

0.1020 

13.582 

Slope  of  hill 
R 

3 

0.3830 

0.1277 

16.991 

AB  interaction 

4 

0.0033 

0.00G8 

0.1090 

AR  interaction 

6 

0.0381 

0.0064 

0.8450 

BR  interaction 

6 

0.0165 

0.0028 

0.3670 

ABR  interaction 

12 

0.0091 

0.0008 

0.1000 

Error 

179 

1.3449 

0.0075 

Each  of  the  three  main  effects  above  are  significant  at  the 
0.01  level,  but  no  interaction  terms  are  significant. 

c.   Constant  and  Dependent  Speed. 

The  ANOVA  for  comparing  constant  and  dependent 
speed  is  given  in  Figure  10.   The  treatment  level  data  is 
as  follows: 

Mean     Standard  Deviation 

Constant  Speed      0.569  0.1732 

Dependent  Speed     0.5393         0.1921 
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2.   Number  of  Intervisibility  Segments 
a.   Constant  Speed 

SUMMARY  TABLE  -  ANOVA 

1/2 

(X    Transgeneration  Data) 


Source 


DF 


SS 


MS 


F-Ratio 


#  of  hill  2 
A 

Height  of  hill  2 
B 

Slope  of  hill  3 
R 

AB  interaction  4 

AR  interaction  6 

ABR  interaction  12 

Error  179 


8.3415  4.1708 

0.2402  0.1201 

2.7278  0.9093 

0.0659  0.0165 

0.2686  0.0448 

0.1389  0.0116 

9.0475  0.5056 


82.281 

2.369 

17.938 

0.3255 
0.8838 
0.2288 


From  the  ANOVA,  the  number  of  hills  and  slope  are  significant 
at  the  0.01  level,  but  the  height  of  hills  is  not  significant, 
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b.   Dependent  Speed 


SUMMARY  TABLE  -  AN OVA 


1/2 
(X    Trans generation  Data) 


Source 

DF 

SS 

MS 

F-Ratio 

#  of  hill 
A 

2 

8.2224 

4.1112 

79.337 

Height  of  hill 
B 

2 

0.2561 

0.1281 

2.4272 

Slope  of  hill 
R 

3 

2.6384 

0.8795 

16.972 

AB  interaction 

4 

0.06446 

0.01611 

0.311 

AR  interaction 

6 

0.2905 

0.0484 

0.934 

ABR  interaction 

12 

0.1289 

0.0107 

0.2065 

Error 

179 

9.2756 

0.0518 

0.2065 

The  results  of  the  dependent  speed  case  are  analogous  to 
those  of  the  constant  speed  case.   That  is,  the  number  of 
hills  and  slope  of  hills  are  significant,  but  the  height  of 
hills  is  not  significant. 

The  mean  and  standard  deviation  of  the  treatment  levels 
are  given  below. 
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SLOPE  MEAN 

CONST      DEP 


100 
200 
300 


STD  DEV 
CONST       DEP 


DOWN 

4.2075 

4.2075 

1.04444 

1.0444 

SIDE 

4.1852 

4.1852 

1.1339 

1.1339 

UP 

3.3148 

3.3333 

1.1947 

1.2131 

BOTH 

4.4259 

4.4259 

1.3818 

1.3818 

#  OF  HILL 

MEAN 

STD 

DEV 

CONST 

DEP 

CONST 

DEP 

6 

4.930 

3.038 

1.1871 

0.7166 

12 

4.139 

4.139 

1.039 

1.0388 

18 

3.042 

4.931 

0.7207 

1.1788 

HEIGHT  OF 

MEAN 

STD 

DEV 

HILL 

CONST 

DEP 

CONST 

DEP 

3.875  3.875 
4.069  4.069 
4.155     4.169 


1.393  1.393 
1.191  1.191 
1.191      1.195 
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V.   CONCLUSION  AND  RECOMMENDATION 

Terrain  is  one  of  the  key  factors  in  the  analysis  of 
land  combat.   Thus,  the  representation  of  terrain  in  a  land 
combat  model  is  a  critical  factor  in  the  model  design.   High 
resolution  combat  models  typically  employ  actual  terrain 
data  to  describe  the  terrain  environment.   This  approach, 
although  realistic  for  the  specific  terrain  being  represented, 
is  very  costly  and  time  consuming  to  implement. 

The  goal  of  the  research  by  Needels  and  this  thesis  is 
to  provide  a  parametric  representation  of  terrain  usable 
either  as  input  to  a  high  resolution  combat  model  or  as  a 
"stand-alone"  model  to  investigate  intervisibility  segment 
length  parameters.   In  particular,  the  current  model 
describes  the  terrain  as  a  function  of  the  number  of  hill 
masses,  the  height  of  the  hills,  the  position  and  spread 
of  the  hills,  and  the  slope. 

A  MBVN  distribution  function  and  a  linear  function 
cutting  the  MBVN  function  are  employed  to  achieve  the 
terrain  representation.   The  terrains  given  in  Figures  12 
through  17  illustrate  the  flexibility  of  the  model.   In 
fact,  with  appropriate  selection  of  the  input  parameters 
(given  in  Figure  9)  any  desired  terrain  may  be  generated. 

The  results  of  the  analysis  described  in  Chapter  IV 
indicate  that  slope  is  highly  significant  in  the  resulting 
values  of  intervisibility  parameters .   Some  of  the  trends 
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noted  from  the  results  given  in  Chapter  IV  are  described 
below. 

1.  An  up  slope  (cut  from  target  side  to  observer  side) 
resulted  in  the  largest  value  of  percent  intervisibility 
time. 

2.  Increases  in  the  number  and  height  of  hills  resulted 
in  decreases  in  the  percent  intervisibility  time. 

3 .  Speed  dependent  on  terrain  slope  resulted  in  a 
slight  reduction  in  percent  intervisibility  time  from  the 
constant  speed  case. 

The  analysis  of  the  hypothetical  data  described  in 
Chapter  IV  serves  to  illustrate  the  potential  utilization 
of  the  model.   The  trends  indicated  from  the  analysis  are 
appropriate  for  the  selected  terrains  and  conditions,  but 
should  not  be  construed  as  generalized  results. 

An  alternate  method  (tilted  method)  for  generating 
unsymmetric  terrain  outlined  below  is  suggested  as  an  area 
for  future  enhancement  of  the  model.   Given  a  MBVN  density 
function  centered  at  Y  in  Figure  11,  define  9  as  the  tilted 
angle,  assuming  p  =  0  .   The  height  of  the  density  function 
is  Z2  at  Y2  and  is  given  by 

1  x-y   ~    y-y  0 

Z2   =   MH*EXP{-  4  [(-—£)   +  (— -*)  ])  (38) 

2  ax        ay 

and 

Z2  rlT 


Y2-Y 


=   tan  q   -  6]  (39) 


4-2 


y 


FIGURE  11   TILTED  METHOD 


From  38  and  39, 


mh*exp{-  |  n-^-V  +  (-^n)  =  (y-y2)  tan(J-e) 


(40) 


For  a  specified  value  of  x  and  y,  the  value  of  y~  can  be 
computed  using  the  golden  section  method.   The  tilted  height, 
ZH,  can  then  be  computed  as  follows: 


y-y. 


ZH   = 


sin  9 


(41) 


Additional  investigation  of  this  method  is  required  before 
implementation  in  the  model. 

In  conclusion,  a  model  has  been  developed  and  exercised 
which  provides  a  potentially  valuable  tool  for  efficient 
analysis  of  various  terrain  effects  in  the  land  combat 
environment.   The  intervisibility  parameters  output  from 
the  model  can  serve  as  critical  inputs  to  land  combat  models 
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for  evaluating  such  things  as  the  contribution  of  mobility, 
agibility,  and  firepower  to  battlefield  survivability. 
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FIGURE  12-2. 


Terrain:  1,  #  of  Hills;  18,  Average  Height; 
100,  Slope;  Down 
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FIGURE  13-2 


Terrain:  2,  #  of  Hills;  18,  Average  Height; 
300,  Slope;  -Up-   k 
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FIGURE  14-2. 


Terrain:  2,  #  of  Hills;  12,  Average  Height; 
100,  Slope;  Up 
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FIGURE  15-2 


Terrain:  3,  #  of  Hills  i  6,  Average  Height; 
200,  Slope;  Side 
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FIGURE  16-2. 


Terrain:  3,,  #  of  Hill?j  6_,  Average  Height; 
300,  Slope;-  Side .^ 
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Terrains  "X;;#^Bf  fidns-;  12,  Average  Height; 
200,  Slope T 'Botk"  56 
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